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Skyrmions, once a hypothesized field-theoretical object believed to describe the nature of ele- 
mentary particles, became common sightings in recent years among several non-centrosymmetric 
metallic ferromagnets. For more practical applications of Skyrmionic matter as carriers of informa- 
tion, thus realizing the prospect of "Skyrmionics" , it is necessary to have the means to create and 
manipulate Skyrmions individually. We show through extensive simulation of the Landau-Lifshitz- 
Gilbert equation that a circulating current imparted to the metallic chiral ferromagnetic system can 
create isolated Skyrmionic spin texture without the aid of external magnetic field. 
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I. INTRODUCTION 

Skyrmions are topological configurations of a vector 
order parameter in field-tlieoretical systems, whose exis- 
tence had been anticipated mathematically by the par- 
ticle physicist Tony Skyrme for half a century^. The 
two-dimensional version of it, sometimes known as the 
baby Skyrmion, was recognized to be the topological so- 
lution of the non-linear sigma model, which serves as a 
model for ferromagnets in two dimensions^. The dis- 
covery of two-dimensional Skyrmions took place in a 
variety of condensed matter systems such as quantum 
Hall ferromagnets'^'*, metallic chiral ferromagnets^"^, 
as well as the ferromagnetic monolayer^ and doped 
antiferromanget^^, and is predicted to be possible in 
two-dimensional liquid crystal systems^^. A series of 
metallic magnets of B20 structure including MnSi^, 
Fei_a;Coa;Si^'^, and FeGe^ is now known to host the 
Skyrmionic magnetic texture. In these magnetic mate- 
rials the formation of Skyrmions is a consequence of the 
competitive interplay between the Dzyaloshinskii-Moriya 
(DM) interaction which tends to induce magnetic spirals, 
and an external magnetic field which tends to bring about 
a ferromagnetic background^^'^'^. Initial theoretical con- 
sideration assumed a layered structure of Skyrmions em- 
bedded in a three-dimensional crystal and found a small 
region for its stability just below the magnetic ordering 
transition temperature^'^^. Later, Monte Carlo simula- 
tion for two-dimensional magnetic model showed a much 
wider region for its existence extending down nearly to 
zero temperature^'^*, as was indeed verified in experi- 
ment on a thin-film Feo.sCoo.sSi^, due to the lack of 
competing magnetic structures in the thin-film geometry. 
Upon increasing the thickness of the film, the Skyrmion 
phase gradually moves to a higher temperature regime*. 

Aside from their fundamental scientific allure, it is of 
importance now to ask if these novel topological objects 
can be manipulated to the extent that they can be used 
as carriers of charge, spin, and of information in gen- 
eral. Applying external magnetic field to a chiral ferro- 



magnet generates a lattice of Skyrmions, or a Skyrmion 
crystal^"*. Spin wave fluctuations can be induced in the 
Skyrmion crystal by applying an additional ac magnetic 
field as was recently studied theoretically^^'^^. Other 
than by magnetic field, the Skyrmions can be manipu- 
lated through the Hund's rule coupling of the spin of 
electrons in the metallic host to the localized moments 
forming the Skyrmionic texture. Recent demonstration 
of the rotation of the Skyrmion crystal axis by the ap- 
plied current*^'** is an example of this kind. In the ab- 
sence of pinning, Skyrmions are expected to move with 
the velocity equal to the drift velocity of the electrons re- 
sponsible for the current^^. Gilbert damping modifies the 
Skyrmion trajectory in such a way to induce Hall-like mo- 
tion of Skyrmion orthogonal to the current direction*^ . 

Although various aspects of current-induced motion 
of Skyrmions have been studied theoretically in a num- 
ber of papers already^®'^^, the possibility of generating 
Skyrmions by the current has not been addressed yet. 
From a technological point of view, it is desirable to have 
the extra means to generate and destroy Skyrmions one 
at a time, instead of having to generate them only in the 
form of a crystal by the external field. Such possibility, 
once realized, will help usher the era of "Skyrmionics" 
- an effort to tailor Skyrmions as carriers of informa- 
tion bits. In this article, using extensive simulation of 
Landau-Lifshitz-Gilbert (LLG) equation, we provide ev- 
idence that circulating spin current can couple to local 
magnetic moments to produce Skyrmion spin texture in 
the latter. Basic theory and LLG equation are introduced 
in Sec. II. Numerical findings are given in Sec. III. Pos- 
sibility of experimental realization will be discussed in 
Sec. IV. 



II. FORMULATION 

We imagine a thin slab of chiral ferromagnet, such that 
all degrees of freedom behave identically along the thick- 
ness direction. The dynamics of the localized moments 



n are governed by the Hamiltonian^'^^ 



H„ 



J 






+ — / d^r n • V X n / d^r B • n. (1) 

a^ J w^ j 

Formulated in the continuum, a expresses the linear di- 
mension of the suitably chosen unit cell, while J and D 
are the ferromagnetic and DM exchange energies, respec- 
tively. The Zeeman field B is typically applied perpen- 
dicular to the planar direction, B ~ Bz. The dynamics of 
the conduction electrons and the coupling of their spins 
to the local moment, on the other hand, are embodied in 
the sd Hamiltonian, 



^- = /^^^*^(l^-^-' 



cr • n MP, * = 



(2) 



Electron operators are given in the spinor form ^ com- 
prising up (c-f-) and down (c^) spin components. Deriva- 
tion of the subsequent equation of motion of n based on 
the total Hamiltonian H = H^ + Hgd is well-documented 
in the literaturc'^'^ and reproduced briefly here. 

In the large Hund's coupling limit {Jh — > oo) the spin 
\E''I'<T\I> of the conduction electron is forced to align with 
the local magnetization direction n while the electron 
band with anti-parallel spins forms a higher-energy con- 
tinuum. The idea is implemented by making the unitary 
transformation, U^m ■ crU = Gz-, where U is given by 



where we also introduced paramagnetic spin current j 
and the spin density p as^" 



J = 



1 



2™ 

p = tp^ip 



(^t[p^]_[p^t]^)^ 



(6) 



Hn 



We focus on the dynamics of the local moments using 
-ffa-j as the total Hamiltonian and ask how an 
externally imposed spin current pattern j influences the 
magnetic dynamics. Already the likelihood of Skyrmion 
induction by the circulating j can be seen in the following 
simple argument. In a steady state, the divergence of 
the current vanishes, V • j = 0, making it possible to re- 
write j as the curl j = V x c. For circulating current of 
cylindrically symmetric form, the vector field c is directed 
in the z-direction. Integrating the energy functional by 
parts, one can re-write a • j ^ — c • V x a and find the 
effective magnetic field (V x a)z couples linearly to the 
current source c. This effective magnetic field is nothing 
other than the Skyrmion density through the relation 



(V X a)z = -n • {d^n x dyu). 



(7) 



A suitable c (j) of sufficient strength will overcome the 
nucleation energy cost and induce a Skyrmion with (V x 
a)z y^ 0. The sign dictates that the creation of Skyrmion, 
Q > 0, is energetically preferred for counter-clockwise 
(CCW) direction of spin current, while anti-Skyrmions, 
Q < 0, are preferred for clockwise (CW) spin current 
direction. Here 



U 



z\ 



cos—, Zi 
2' 



e*^ sin ■ 



(3) 



The spinor z = 



forms the CP^ representation of the 



classical spin n — ■z*a'L = [sin 6 cos cf), sm 9 srcKJ)^ cos 0). 
The new electron spinor [/^^ has the upper (lower) com- 
ponent parallel (anti-parallel) to n, out of which we 
choose to keep the upper one only, denoted %!). The field 
operator ip gives the electron whose spin direction is pro- 
jected strictly parallel to the local moment n. The Hamil- 
tonian for f/j is 



H^„ 



d^r V^- 



haf 



2m 



4. 



(4) 



where the gauge field a is derived from z as a^ = —i'z,'^ di'z,. 
It is gauge-invariant under z — >■ e*-^z, and ip — > e~^^ip. 
Expanding the sd Hamiltonian gives the coupling of local 
moment to the electrons as^*^ 



H, 



_,^/.^r(na.j + |^a^ 



(5) 



Q^— d^rn- (d^n x ^j^n) (8) 

is the total Skyrmion charge. 

The real-time action for the magnetic dynamics is 



-^ I d^vdt (1 - cose)dt<p - / dt{Hn + i7a-j) (9) 

with the effective moment 5, and (0, (j)) are the polar and 
azimuthal angles of the magnetization unit vector n. The 
LLG equation follows as 



h + -n X (- Ja^V^n + 2DaV x n - B) 
n 

+ an X h + — (j • V)n = 0. 



(10) 



The Gilbert damping constant a is introduced 
phenomenologically^". Assuming all the spins sharing 
the same (x, y)-coordinate behave identically, we arrive 
at a two-dimensional discretized version of the LLG equa- 
tion 



hj + Ilj X 2_](-ni+aej. + nrXi+aej, X Cji - K^h 



jei 



an," X n,- 



V^iUi+ax - Hi) + wf (rii+ay - Hj) ) = 0. (11 



The drift velocity v^ — (ff ,wf) is related to the current 
density through a^j — pxw, where x is the number of con- 
duction electrons enclosed in a unit cell of volume o'^ and 
p is the spin polarization fraction. All physical quantities 
in the equation are made dimensionless by choosing the 
unit of time to = U/J, and k = D/J, b = B/{D'^/J). 
The quantity vq = aS/topx with the dimension of veloc- 
ity is introduced as well. The wavelength of the spiral 
is given by A ^ 2n\/2a/K in two dimensions. Square 
lattice of size L x L is adopted throughout the simula- 
tion and the k value equal to 0.5, which yields the spiral 
wavelength A ^ ATT^/2a ^ 17.8a. The diameter of the 
Skyrmion is roughly A/2, which makes the Skyrmion ra- 
dius Rs ^ A/4. The sum X^isi spans the four nearest 
neighbors of the site i, connected by the unit vectors Cji. 
Time- integration of the discrete LLG equation (11) is 
done using the fourth-order Runge-Kutta method. Each 
unit of time to is sliced into over 1000 steps for the 
integration procedure, with no difference in the result 



J 



found for finer step sizes. Typical integration ran up to 
t/to ~ 10"^ without a noticeable error accumulation. The 
spin current density j = v/wq is given the circulating 
profile with the Lorentzian shape 



j(r) = 4>jo 



R,. 



i?? 



(12) 



at the center of the simulation lattice, characterized by 
its magnitude jo and the extent Re- Both positive and 
negative jo corresponding to the CCW (CW) circulation 
of the spin current were considered. In our convention, 
a uniform spin current j = jo leads to the Skyrmion flow 
velocity Vs in the same direction, v^ || jo, opposite to 
what was used in Rcf. 19. Previous study^'^^'^^ found 
the magnetic phase with helical spin configuration at zero 
magnetic field and Skyrmion crystal phase over the inter- 
mediate field values bd ^ b ^ bc2- Two magnetic fields 
are considered in our simulation here, one for b slightly 
less than the upper critical field strength 6c2 where a typ- 
ical ground state consists of ferromagnetic spins dotted 
with a few isolated anti-Skyrmions (Fig. Ib)^'^'^*^, and 
the other at 6 = and having the helical spins as the 
ground state (Fig. 2b). We refer to these two ground 
states as FM' (FM-|-anti-Skyrmions) and H, respectively. 



r 



III. PROCESSES OF CURRENT-INDUCED 
SKYRMION FORMATION 

A. Skyrmion production 

Due to the simpleness of the ferromagnetic spin back- 
ground, it is easier to understand the formation process 
of Skyrmions for the FM' case, b < bc2- For optimal 
Skyrmion generation condition we choose Re comparable 
to the typical radius of a Skyrmion, Re « Rs- As the 
system evolves over time, the time-dependent Skyrmion 
number Q{t) = (l/Air) J (fr n(t) • {d^n{t) x dyn{t)) is 
evaluated to keep track of the Skyrmion creation process. 

As shown in Fig. 1, CW spin current produces, for a 
brief moment, a spin texture that corresponds roughly to 
a Skyrmion-anti-Skyrmion pair within the radius r '^ Rs 
(red core region in Fig. Ic-d). Out of the pair, the 
Skyrmion part soon vanishes, leaving behind an anti- 
Skyrmion (Fig. le-g). The anti-Skyrmion then gradually 



drifts out, due to the finite Gilbert damping which gives 
rise to a velocity orthogonal to the local current flow. 
The Gilbert damping-induced Skyrmion Hall motion was 
extensively discussed in Ref. 19 to which we refer the in- 
terested readers. The Skyrmion meanwhile experiences 
an inward flow and vanishes. Once the anti-Skyrmion 
has drifted a sufficient distance out, another burst of 
Skyrmion-anti-Skyrmion pair takes place (Fig. Ih), with 
the Skyrmion portion vanishing through the core again 
decreasing the total Skyrmion number by an integer 
amount (Fig. li-1). Each burst of the Skyrmion-anti- 
Skyrmion pair and the subsequent decay of the Skyrmion 
is responsible for the total Skyrmion number Q{t) jump- 
ing by an integer amount as shown in Fig. la. 

Evidently creation of Skyrmion-anti-Skyrmion pair 
(Skyrmion dipole) is energetically cheaper than creat- 
ing an isolated Skyrmion (Skyrmion monopole). Once 
a pair is created, the Gilbert damping-induced Hall mo- 
tion naturally moves the two objects in the opposite di- 
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FIG. 1: (color online) Skyrmion generation by circulating spin current source in the FM' background, (a) Time dependence 
of the total Skyrmion number Q{t). Q(0) ~ —2 is nonzero from the residual anti-Skyrmions in the ferromagnetic background. 
Blue (red)-colored curves correspond to CW (CCW) circulating current with the jo value indicated for each curve, (b)-(l) 
Time-dependent snapshots of the spin configuration over the green-circled time interval in (a). Adjacent spins are separated by 
the distance a. The circle at the center of (b) indicates Re, the extent of the current source. Residual anti-Skyrmions lie outside 
the field of view. Red and blue refer to positive and negative Skyrmion densities, respectively. Time at which the snapshot is 
taken and the change in the Skyrmion number are given as (i/io, Q(i) — Q(0)) above each figure. Skyrmion number is seen to 
decrease by one, and later by three. Threshold value joc lies between 2.35 and 2.4 for CW spin current. 



rections, so that only one species of Skyrmions can leave 
the current core. When CCW current (jo > 0) is ap- 
plied, the same Gilbert drifting mechanism that pushed 
the anti-Skyrmions out now tends to absorb the nearby 
anti-Skyrmions toward the core region because the di- 
rection of the Hall motion is also reversed. When an 
anti-Skyrmion is pulled sufficiently close to the center, 
a brief burst of Skyrmion appears at the core and "eats 
up" the attracted anti-Skyrmion. A good mental picture 
of the process is to take Fig. Ib-lg and reverse them 
in time. When all the existing anti-Skyrmions are eaten 
up in this way, the total Skyrmion number Q{t) remains 
close to zero (Fig. la) and does not rise above it as a 
consequence of the fact that positive charge Skyrmions 
are energetically forbidden for the magnetic field 6 > 0. 
There is a critical value of the spin current density 
joc required to produce anti-Skyrmions out of the ferro- 
magnetic spin background. Obviously the rate of energy 
input from the circulating current needs to outpace that 
of the energy drain due to the Gilbert damping in order 
to supply sufficient energy to create an anti-Skyrmion. 
Above the threshold, |jo| > joc, the total Skyrmion num- 
ber begins to jump in integer steps with time intervals 
that decrease as |jo| — joc increases (Fig. la). 



Figure 2 depicts the Skyrmion creation process for H 
background (6 = 0). As with the FM' background, Q{t) 
can be seen to increase in integer quantized steps. In 
contrast to & > case which prefers an anti-Skyrmion 
charge, opposite signs of the circular current result in 
more or less symmetrical profiles of Q{t) as shown in Fig. 
2a. Snapshots of spin configurations for Q = — >■ Q = 1 
are in Fig. 2b-g, and those ofQ^l— >Q = 2in Fig. 2h- 
21. At first a small patch of helical spin is torn into two 
disjoint segments as in Fig. 2c-d near the current core, 
giving way to the nucleation of some defects in the in- 
tervening region. The defect profile is best described as 
the meron-anti-Skyrmion-meron composite (merons are 
in blue, ant-Skyrmion in red in Fig. 2c), out of which 
the central anti-Skyrmion shrinks in size while the two 
merons drift out, again due to the Gilbert damping- 
induced Hall motion. The two merons are responsible 
for the total charge Q{t) ss 1. Each time the process 
repeats itself, the Skyrmion number increases by one. 

In both H and FM' spin backgrounds, the resulting 
charge of the Skyrmions produced by the current is cor- 
related with the helicity of the circulation as predicted 
in an earlier argument; see paragraph following Eq. (6). 
In detail, however, the production process proves to be 




b.(l,0) 



c. (50, -0.01) 



d. (130, 0.08) 




e. (142, 0.30) 






,0 92) 


g (155, 93) 


#1*" 


m J^ 



250 

h (500 89) 



500 



750 



1000 




I (680 90) 



J (795, 1 23) 



k (799 1 89) 



Sfe|f «fej|- 






ii 



I. (815, 1.86) 








FIG. 2: (color online) Skyrmion generation by circulating spin current source in the H background, (a) Time dependence of the 
total Skyrmion number Q(t) for various jo values. Initial Skyrmion number is zero for the helical spin. Same color convention 
as in Fig. 1 applies throughout this figure, (b)-(l) Two sets of snapshots corresponding to Q = — s- 1 and Q = 1 — s- 2 
corresponding to two green circles in Fig. 2a. 



far more intricate involving an Skyrmion-anti-Skyrmion 
pair at the initial stage of metamorphosis rather than 
the simple creation of an isolated Skyrmion out of the 
vacuum. A different mechanism then intervenes, namely 
Gilbert damping-assisted Hall motion, which separates 
the pair according to their respective charges, leaving 
only one species of Skyrmions to survive at the end. 



B. Influence of Gilbert damping 

Gilbert damping plays a major role in the generation 
of Skyrmions and ant i- Skyrmions by facilitating their ra- 
dially outward motion from the circulating spin current 
region. All results shown in Figures 1 and 2 are obtained 
for a = 1, which is a fairly large value for Gilbert damp- 
ing parameter. As the damping is reduced, the sharp, 
integer quantization of Skyrmion number Q{t) gets less 
conspicuous as shown in Fig. 3. According to our numeri- 
cal observation, clear integer quantization oi Q{t) is more 
or less synonymous with the production of well-isolated 
Skyrmions in real space, which only becomes possible for 
a in excess of certain minimum value as shown in Fig. 
3b and 3d. This is related to the quite complex manner 
in which the Skyrmion creation process takes place. The 
initial creation process, as shown in Figs. 1 and 2, always 



involves a pair of Skyrmions of opposite charges, out of 
which only one is pushed outwards by virtue of Gilbert 
damping- assisted Hall motion. A large Hall motion re- 
quires a correspondingly large Gilbert damping constant. 
Otherwise, when a is small, those Skyrmions after gener- 
ation continue to move in circular path (because current 
is circular) rather than drifting out, and collide/merge 
with other Skyrmions. This leads to a large overlap be- 
tween spatially adjacent Skyrmions, as well as a non- 
integral value of Q{t)'^'^. Only with sufficiently large a 
do we find that the Skyrmions are being pushed out fast 
enough after their generation to minimize the overlap 
and give rise to integer Q{t). Mochizuki's recent LLG 
simulation^® uses the Gilbert constant a < 0.04, while 
extra damping effects due to spin-motive force are ex- 
pected to increase this value to around 0.1^^. Further 
enhancing the Gilbert damping might require intentional 
disordering of the film by irradiation. 



C. Skyrmion lifetime 

Once the Skyrmions have been generated by means of 
current, it is important that they remain there after the 
current is turned off if they are to serve as memory bits. 
After the current was applied long enough to nucleate a 
Skyrmion, we turned it off in the simulation to see if it 
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FIG. 3: (color online) Influence of Gilbert damping on 
Skyrmion generation, (a)-(b) Plot of Skyrmion number Q{t) 
with & = (H background) for different Gilbert damping 
constants a. |Q(t)l increases monotonically with a without 
showing signs of quantization when a is small. At a larger 
damping, 0.3 < a < 0.5, quantization of Skyrmion number 
begins to occur. The time span of each integer plateau be- 
comes longer, and the production rate of Skyrmions becomes 
less, with increasing a. (c)-(d) Plot of the Skyrmion num- 
ber Q{t) with FM' background. Numbers in the inset are the 
Gilbert damping constants, {jo, Re) = ( — 1.5,3.0) was used 
for (a)-(b) and {jo, Re) = (-2.5,3.0) for (c)-(d). 
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FIG. 4: (color online) Skyrmion generation by pulse cur- 
rent. (a)-(b) Pulse current of size jo = 5 is applied for 
< ti/to < 50 over the helical spin state and then turned 
off. (c)-(d) jo = 15 for < ti/fo < 15 over the FM' state. 
Dotted vertical line in (b) and (d) indicates the termination 
of the applied current. Top row: the spin configuration and 
Skyrmion number Q{t) — Q{0) when a sufficiently long time 
has elapsed after the current is turned off. Bottom row: time 
dependence of the Skyrmion number Q{t) — Q(0). All current 
directions are CCW with Re = 3 and a = 0.1. 



would decay. In fact, the decay of Skyrmion is expected 
because without the current, the Skyrmion would be in 
a metastable state. As the Q{t) plots in Fig. 4b and 4d 
show, however, the ant i- Skyrmions once created remain 
extremely stable although energetically it cannot exist 
in the ground state for 6 = 0. As far as our numerical 
simulation lasts, the spin patterns shown in Fig. 4a and 
4c remained persistently the same. Such a long lifetime of 
an isolated Skyrmion is a positive feature in utilizing the 
chiral ferromagnet as a platform for encoding Skyrmion- 
based information. 



D. Effects of radial current pulse 

Earlier we used the circulating current as a platform 
to generate Skyrmions electrically. Clear, integer quan- 
tization of Skyrmion charge was possible only for very 
large Gilbert damping parameters. Here we show that 
Skyrmions can also be created by radial current pulses 
as given by 



j(r) 



; Jo 



(13) 



While the radial current source remains on, a bunch of 
plus and minus Skyrmion fragments appear from the cen- 
ter and drift out as they collide and annihilate with oppo- 
sitely charged fragments as shown in Fig. 5a. When the 
radial current is turned off, most of the plus and minus 
Skyrmion fragments annihilate, leaving behind a portion 
that evolves into an integer-valued Skyrmion, with the 
sign as favored by the direction of the applied magnetic 
field, as shown in Fig. 5b and 5c. This way of creating 
Skyrmion is only possible when external magnetic field is 
acting on the helical magnet. For helical spin configura- 
tion, the plus and minus Skyrmion densities created by 
the radial current annihilate and disappear completely, 
returning back to the initial helical spin configuration 
over time. 

We emphasize that this method of generating 
Skyrmions in the FM' background works very well even 
for small values of the Gilbert damping constant a. This 
is because the selection process to isolate one particular 
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FIG. 5: (color online) Skyrmion creation by radial current 
source. Current pulse with jo = 10 lasting for t/to — 10 was 
given. The numbers above each figure indicate {t/to,Q{t) — 
Q(0)). (a) Right at the turn-off, one finds a swarm of posi- 
tive and negative Skyrmion fragments, (b) After some time 
has passed, most of the Skyrmions and anti-Skyrmions have 
annihilated each other, leaving behind a predominant anti- 
Skyrmion density due to the external magnetic field, (c) Final 
state shows one anti-Skyrmion. 



sign of Skyrmion is done by external magnetic field, not 
by Gilbert damping as in the case of circulating current. 



IV. EXPERIMENTAL REALIZATION AND 
DISCUSSION 

The typical exchange energy scale of a chiral ferromag- 
net is Jo ~ 3 meV^"^, with the dispersion huk ^ Jo(A;ao)^ 
involving oq, the linear dimension of the physical unit 
cell. On the other hand, from our continuum LLG equa- 
tion in (10) follows the dispersion relation TiuJk ^ J{ka)^. 
Hence, J used in our model is related to the microscopic 
parameter Jq hy J ^ Jo(ao/a)^, which also sets the 
time and the velocity scales as to = Ti/J w 220(a/ao)^ 
fs and vq — aS/topx ^ (ao Jo/?i)(S'/p)(ao/a)'* ^ 1.8 x 
10^{S/p){ao/a)^ m/s. 

In order to fix the ratio a/ao, we note that with our 
choice K = 0.5 the typical wavelength of the spiral is 
A = 4:TTy/2a ss 17.8a. Values of A vary quite markedly 
among the chiral ferromagnets, ranging from ~3 nm for 



-20 nm for MnSi^ ~70 nm for FeGe^ and -90 



FeGe-' 

nm for Fei_3;C0a;Si'''^. Taking the characteristic spi- 
ral wavelength A ~ 60 nm yields 17.8a ~ 60 nm, or 
a ~ 4 nm. Choosing the microscopic lattice constant 
ao — 4A then gives the ratio a/ao — 10. Finally, with 



the spin polarization fraction p '^ O.l and S — 1 we ob- 
tain uq ~ 1 m/s. The critical current density joc of order 
unity we saw in the simulation thus corresponds to the 
drift velocity of a few m/s to be imparted to the con- 
duction electrons. Estimate of the carrier concentration 
in MnSi^* pe -- 10^^ - lO^^/m^ yields the critical cur- 
rent density jc ~ 10^° — lO^^A/m^, which is close to, 
or below the critical current required for domain wall 
switching in typical magnetic devices^'^. The time inter- 
val between successive Skyrmion creation events as shown 
in Figs. 1 and 2 is ^ lOOtg. With the estimate to ^ 20 
ps, the time interval is ~ 2 ns. The requirement for 
Skyrmion nucleation is reduced to having a strong cur- 
rent pulse of a few nanosecond duration with the density 
jc ^ 10"'^° — lO-'^^A/m^ - a practice readily available in 
modern-day spintronics laboratories. 

Although a detailed practical scheme to implement the 
circuitry needed for the Skyrmion generation is yet to be 
worked out, our simulation provides a proof-of-concept 
demonstration for Skyrmion creation by electrical cur- 
rent. It is conceivable, in principle, that an artificial 
bend in the current pathway (like a curved freeway) is 
fabricated to mimic the circular pattern we assumed in 
the simulation. As the current of enough intensity passes 
through such a bend, Skyrmions will "pop out" from the 
point of greatest curvature. The Skyrmion production 
by means of radial current discussed in Sec. Ill D can 
be mimicked by designing a narrow channel of metallic 
magnet that opens up into a much wider region. A cur- 
rent coming through such constriction will be a source of 
radial motion, thus of Skyrmions. Although much engi- 
neering ideas need to be carved out yet, the fundamental 
processes revealed by the present numerical simulation 
will be exciting to realize in laboratories in the near fu- 
ture. 
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